Relevance of multiband Jahn- Teller effects on the electron-phonon interaction in A^C 
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Assessing the effective relevance of multiband effects in the fuUerides is of fundamental importance 
to understand the complex superconducting and transport properties of these compounds. In this 
paper we investigate in particular the role of the multiband effects on the electron-phonon (el-ph) 
properties of the t\u bands coupled with the Jahn- Teller intra-molecular Hg vibrational modes in 
the Ceo compounds. We show that, assuming perfect degeneracy of the electronic bands, vertex 
diagrams arising from the breakdown of the adiabatic hypothesis, are one order of magnitude smaller 
than the non-crossing terms usually retained in the Migdal-Eliashberg (ME) theory. These results 
permit to understand the robustness on ME theory found by numerical calculations. The effects of 
the non degeneracy of the iin in realistic systems are also analyzed. Using a tight-binding model 
we show that the el-ph interaction is mainly dominated by interband scattering within a single 
electronic band. Our results question the reliability of a degenerate band modeling and show the 
importance of these combined effects in the AzCwi family. 

PACS numbers: 74.70.Wz,63.20.Kr, 74.25.Kc 



Although there is a large consensus about the phonon- 
niediated nature of the superconducting coupling in ful- 
lerides AsCgo compounds, the precise identification of the 
origin of so high critical temperatures {Tc up to ~ 40 K) 
is still object of debate. Electron and phonon properties 
in these compounds are described in terms of molecu- 
lar crystals, where the small hopping integral t between 
nearest neighbor Ceo molecular orbitals sets the scale 
of the electronic kinetic energy, while electron-phonon 
coupling is essentially dominated by the intra-molecular 
vibrational modes l]. First-principle calculations sug- 
gest that el-ph coupling is mainly dominated by the Hg 
modes, having energies loh '--^ 30 — 200 meV, with some 
small contribution from the Ag modes. A simple appli- 
cation of the conventional Migdal-Eliashberg theory in 
these materials is however questioned for different rea- 
sons. On one hand, the small value of t gives rise to 
strong electronic correlation effects since t appears to be 
of the same order of the intra-molecular Hubbard repul- 
sion U flj. On the other hand, nonadiabatic effects are 
also expected to be relevant since the electronic energy 
scale t is of the same order as well of the phonon fre- 
quency scale Wph |3|. In addition, the multiband nature 
of the electronic states iiu involved in the pairing and 
the Jahn- Teller nature of the Hg phonon modes make 
the problem even more complex but also more interest- 
ing. 

A popular tool to investigate el-ph properties in sys- 
tems where all the energy scales t, U, Wph are of similar 
order of magnitude, is the Dynamical Mean Field Theory 
(DMFT) 3] which does not rely on any small parameter 
expansion. In DMFT the lattice Hamiltonian is mapped 
in a local impurity model where the electron dynamics 
through the crystal is taken into account by a set (in 
multiband systems) of frequency-dependent Weiss fields. 



This approach is particularly suitable in fuUerides be- 
cause the intra-molecular el-ph interaction, just as the 
Hubbard repulsion, is local by construction when ex- 
pressed in the basis of molecular orbitals. In order to pro- 
vide a close set of equations, the hopping integrals t™™ 
are often chosen to be diagonal with respect to the in- 



dexes of molecular orbitals ' 



and, in the absence of 



crystal field spHtting, degenerate i^™ = Smm'Uj- By us- 
ing this approximation, the interplay between electronic 
correlation and el-ph coupling mediated by Jahn- Teller 
mode has been studied, and a reduction of the effective 
superconducting pairing for small U |^ accompanied by a 
remarkable enhancement of it close to the metal-insulator 
5] transition has been found. Within the same approxi- 
mation the validity of ME theory for a multiband Jahn- 
Teller {t x H) model has also been investigated, showing 
a good agreement between ME theory and DMFT re- 
sults up to intermediate-large values of the el-ph coupling 
A ^ 1, in contrast to the case of a single band system in- 
teracting with a non Jahn- Teller mode (a x A) where ME 
theory breaks down for A ^ 0.5 [J|. The numerical solu- 
tion of DMFT equations makes hard however to identify 
the physical origin of these results. 

Aim of this paper is to provide an analytical insight 
about the relevance of the multiband and Jahn- Teller- like 
effects on the el-ph interaction in fuUerides. In particu- 
lar we show that, assuming degenerate electronic bands, 
the peculiar symmetry of the Jahn- Teller Hg modes of 
the Ceo molecules gives rise to a drastic suppression of 
the el-ph scattering channels described by vertex dia- 
grams, enforcing the non-crossing approximation (NCA) 
which is at the basis of ME theory. Such a suppression 
is even more strict in a two band e x E model where 
vertex diagrams exactly cancel out. At the same time, 
we show that these results strongly depend on the de- 
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FIG. 1: Skeleton diagrams for: (a) el-ph self-energy in non- 
crossing approximation (ME); (b) lowest order expansion of 
the el-ph vertex function. Solid lines are electron Green's 
functions, wavy lines are phonon propagators, and small filled 
circles the el-ph Jahn- Teller matrix elements. Each pictorial 
element is a matrix in the space of the orbital molecular index. 
The same diagrammatic expressions hold true even in the 
superconducting state in terms of Nambu representation. 



generate band assumption, whereas a careful analysis of 
the electronic structure of the A^Ceo compounds points 
out that the effective non-degeneracy of the ii„ bands 
significantly reduces the multiband Jahn- Teller effects. 

The simplest vifay to consider multiband Jahn- Teller 
effects in fullerides is to consider electrons in a tight- 
binding model on a fee lattice interacting with a single 
five-degenerate Hg phonon mode with frequency ujq. The 
Hamiltonian reads thus: 

m^m' — l ij M— 1 ^ 

5 3 

^ M— 1 ra,m' — l i 

where m, m' are indexes of the molecular orbitals, i is the 
site index and /i labels the five-degenerate vibrational 
modes, c (c^) and a (a^) are the usual creation (an- 
nihilation) operators for electron and phonon with the 
quantum number specified in the Hamiltonian, and the 
matrices Vj^^, are reported in literature d| . This model 
will be compared in the following with a single band non 
Jahn- Teller model defined as: 
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The prefactor I/vTO in the el-ph scattering term of Eq. 
(jT]) has been introduced in order to normalize the effective 
coupling with Eq. |5J, namely: 
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Let us consider for the moment the case of diagonal 
and degenerate hopping integrals t™™ — 5mm' tij^ as it is 
usually done in DMFT. In this case the molecular orbital 
index m represents a good quantum number also for the 



Bloch-like electron bands, which remain degenerate, and 
both electron Green's function G and its self-energy E 
act thus as the identity matrix / in the molecular orbital 
index G — GI, E = E/. The el-ph matrix interaction at 
the ME level (Fig. ^) can be evaluated in the m x m' 
(3 X 3) multiband space, and it reads thus: 



^ME 




(4) 



where Xq = 2g'^N{Q)/oJo is the single degenerate contri- 
bution from intraband and interband scattering and A^(0) 
is the electron density of states (DOS) for a single band 
of the three-degenerate model) . The effective el-ph cou- 
pling in the normal state self-energy and in the Cooper 
channels are respectively given by X^ = J2m' ^mm' = Aq 
and by the maximum eigenvalue of the matrix Q Age — 
maxcig[A] = Ao |3, just as in the simple single band anal- 
ysis. As a matter of fact, it is easy to see that the matrix 
Q can take a block diagonal form and it behaves as a 
non Jahn- Teller single band system with effective el-ph 
coupling Ao in the one-dimensional symmetrical subspace 
c = l/\/3 X^m cJralO)- Note however that some residual el- 
ph scattering with coupling Aq/I is still operative in the 
orthogonal non-symmetric space in the molecular orbital 
index, although it has no contribution at the ME level 
if the band degeneracy is preserved. As we are going to 
see, this small component with el-ph coupling Aq/IO is 
the only contribution surviving in the higher order terms 
which involve vertex diagrams beyond ME theory. 

We discuss now the robustness of the ME theory with 
respect to the inclusion of highest order diagrams, i.e. 
vertex corrections, which are not taken into account in 
the ME framework. These diagrams are usually neglected 
in the ME theory for single band systems by virtue of 
the so-called Migdal's theorem which states that el-ph 
interaction processes containing vertex diagrams, as in 
Fig. ^, scale with the adiabatic parameter cuo/t and are 
thus negligible in conventional materials where too/t <C 1 
[SJ. This results does not apply a priori however in the 
case of Cgo compounds and of other narrow band sys- 
tems where ljq ^ i |9j, so that the agreement between 
numerical DMFT data and ME theory remains somehow 
surprising. As we are going to show, the physical ori- 
gin of such agreement stems from the particular matrix 
structure of the V matrices in the Jahn- Teller multiband 
case of the fullerides which leads to an almost complete 
cancellation of the vertex diagrams independently of the 
value of the adiabatic parameter LOo/t. 

Once more, the simplest way to understand this fea- 
ture is to assume a diagonal-degenerate hopping term. In 
order to clarify the role of multiband Jahn- Teller effects 
on the vertex processes, let us compare at the skeleton 
level, for both t x H and a x A models, a self-energy 
contribution involving a vertex diagram (Fig. ^p) and 



a typical vertex-free diagram which is usually taken into 
account in ME theory (Fig. ^). For a, a x A model we 
can write in a compact form: 

S°(fc) = -ff^^i?(fc-p)G(p), (5) 

p 



where p is the momentum-frequency vector (p,a;p) and 
^ = '^r, J dujk/'^'^, and where A(p, A:) represents the 
lowest order vertex correction depicted in Fig. ^. In the 
early 60's Migdal was able to show that 



lim A{p,k)(x—, 



(7) 



which allows to neglect vertex diagrams in the adiabatic 
limit uJo/t <^. 

The same diagrammatic picture (Fig. ^ and similar 
analytical expressions hold true for the multiband Jahn- 
Teller t x H model, properly generalized in the matricial 
space of the molecular orbital index. For the diagonal- 
degenerate hopping term case the electron Green's func- 
tion are once more proportional to the identity matrix 
/, so that the only matricial structure comes from the V 
matrices in the el-ph scattering term. In explicit way we 
can write: 



sr(fc) = -^ 
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where the last relation comes from the matricial proper- 
ties of the V terms: E^^ V''V''V''V'']/100 = I/IO. 

The comparison between Eqs. (0-® and Eqs. Q-® 
shows that a strong reduction of the vertex diagrams, of 
a factor 10, is operative in the t x H model appropri- 
ate for fuUerides compounds in the case of a diagonal- 
degenerate hopping term, validating at a large extent 
the non-crossing approximation well beyond the adia- 
batic regime. We would like to stress that the robustness 
of the ME theory with respect to the inclusion of ver- 
tex diagrams does not stem from the negligibility of the 
vertex function A, which can be well sizable for luq ~ i, 
but from the non commutativity of the V matrices and 
from the assumption of diagonal and degenerate hopping 
term. As matter of fact, these results hold true even if 



the assumption of a diagonal hopping term is relaxed as 
long as the three electronic bands are degenerate. To 
show this, we can diagonalize the hopping term in Eq. 
(|T|) by the transformation: 



Ca(k) 



j2Mrc 
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SO that the el-ph scattering term becomes: 



(10) 



(11) 



Since the electron Green's function are still diagonal and 
degenerate G{k) — IG(k) in this new basis, the matricial 
structure of the self-energy terms is still only given by 
the el-ph matrices U. We have thus 
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Finally, it is easy to check that, if the double-degenerate 
Jahn- Teller e x E model is considered, the non- 
commutativity of the corresponding V matrices leads to 
a complete cancellation of the vertex diagrams. 

From the above discussion, one could be tempted to 
conclude that an effective non-crossing ME theory is ex- 
pected to be enforced in fuUerides due to an almost com- 
plete cancellation of the vertex processes. However, it 
should be stressed the above results are valid as long 
as the electronic bands can be considered degenerate. 
In the last part of this paper, we are going to argue 
that this latter assumption is not appropriate in A3 Geo 
compounds in regards with el-ph properties, and a sin- 
gle band model, which is unaffected by Jahn- Teller ef- 
fects, is a better starting point for uiq <C t. In order 
to introduce a more realistic electronic band structure 
than a simple diagonal-degenerate model in the molecu- 
lar orbital space, we consider a tight-binding (TB) model 
for the tiu bands as discussed in Ref. [lOJ which repro- 
duces first-principle LDA calculations with a high degree 
of precision. For sake of simplicty we consider the one- 
directional TB model where all the Cgo molecules have 
a fixed orientation in the fee cystal. Different directional 
ordering will be discussed in a more extended publica- 
tion but they are not expected to qualitatively affect our 
results since they still predict a set of narrow bands with 
only few Fermi cuts. The electronic hopping term in Eq. 
(^ is then diagonalized by the M matrix introduced in 
(|10|l . where each matrix element depends on the elec- 
tronic momentum k. For a generic point of view, the de- 
generacy of the electronic states in the crystal structure 
is preserved only on special high-symmetry points. The 
electronic band structure along the high-symmetry axes 
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FIG. 2: (color online), (a) Total DOS (black solid line) and 
partial DOS for each electronic band in the one-directional 
TB model for A3 Ceo- The dashed vertical line represents 
the Fermi level /i. (b) Left panel: Ratio between the only 
intraband contribution A22 and the total el-ph coupling Ame. 
Right panel: band components \cif of the eigenvector of the 
multiband el-ph coupling matrix A. 



of the fee Brillouin zone and the total electron DOS has 
been widely reported in literature [1|. In order to assess 
to which extent the realistic band structure can be ap- 
proximated by a degenerate three-band model, it should 
be noted that (low energy) electronic and transport prop- 
erties are mainly determined by quantities defined at the 
Fermi level, as the Fermi surface itself, the electron DOS 
at the Fermi level, the Fermi velocity, etc. A realis- 
tic band structure can thus be modeled as a degenerate 
three-band system only if it presents qualitatively similar 
Fermi sheets, with similar densities of states and similar 
Fermi velocities. In this perspective A^Cqo compounds 
seems not to be a good candidate since it is known to 
have only two Fermi cuts of the electronic bands with 
quite different Fermi surfaces [Ij . The smooth topolog- 
ical evolution of the Fermi surfaces with the Fermi level 
permits to identify in an unambiguous way, for what con- 
cerns electronic and transport properties, the three elec- 
tronic bands with their corresponding DOS. The elec- 
tronic DOS for each band, as well as the total DOS, is 
plotted in Fig. |2^, which shows that the electronic band 
structure is largely dominated at the Fermi level // = — 23 
meV, by a single band which is roughly half-filled, with 
two additional bands which are quite far from the Fermi 
level. Only two bands have finite DOS at the Fermi level, 
in agreement with the report of only two Fermi surface. 



We would like to stress that the non-degeneracy of the 
DOS is not related to the removing of the crystal field 
symmetry, but it is uniquely determined by non degen- 
eracy of the electronic dispersion. If we evaluated the 
mean value e = (e) and the variance Se = ((e — e)^)^/^ 
of each band, we obtain ei — —120 meV, Sei = 60 mcV, 
£2 = -23 meV, Se2 = 47 meV, €3 = 135 meV, ^£3 = 30 
meV. This considerations suggest ^3 Ceo can be modeled 
in good approximation as a single band system and that 
multiband effects are qualitatively irrelevant as far as the 
phonon energy is not sufficiently high to switch on el-ph 
scattering between different bands. 

We sustain this intuitive argumentations with a nu- 
merical calculation of the multiband el-ph matrix inter- 
action Ak within the non-crossing approximation (Fig. 
^) which generalizes the ME theory to systems with 
generic DOS N{e). For sake of simphcity we consider 
unrenormalized electron and phonon propagators, with 
energy spectra respectively given by the TB electron 
dispersion and by the phonon frequency loq. A k- 
independent el-ph interaction is obtained in the usual 
way by a proper average over the Fermi surface: A"'' — 
Ek^K)Af]/Ek^(ek)], where D(e) = (lA)c.o/K' + 
e^), which reduces to J2k ^i^k) — -^(^f) ^O'" ^0 -^ 0. For 
sake of simplicity we neglect self-energy effects, which are 
expected to shrink furthermore the electron bands at the 
Fermi level. In order to provide a first simple estimate 
of multiband effects, in the left panel of Fig. ^jp we com- 
pare the effective total el-ph coupling Ame, obtained as 
the maximum eigenvalue of the matrix A, with the single 
intraband contribution A22 relative to the high DOS cen- 
tral band shown in Fig. ^. For ujq -^ 0, A22/AME — 0.86, 
pointing out that the dominant el-ph scattering comes 
from intraband processes with a small interband contri- 
bution due to the small DOS of the lowest electron band. 
As expected the discrepancy between A22 and Ame be- 
comes more relevant as the activation of high energy in- 
terband processes by the phonon increases. To quantify 
further the role of multiband effects we plot in the right 
panel of Fig. |3) the components of the eigenstate of the 
el-ph scattering matrix on the Bloch-like band index a. 
For ujQ ^ only the |c2p component is sizable reflecting 
that cl-ph scattering is mainly intraband. By increas- 
ing Wo interband processes are gradually turned on, but, 
for the physical range of wq ^ 200 meV, intraband el-ph 
scattering within the central band are still dominant. In 
the opposite limit uj ^ W (W being tht total electronic 
bandwidth) the three components become identical and 
the A^Ceo compounds are expected to behave has a de- 
generate three band system. 

In conclusion, in this paper we have investigated the 
role of the electronic structure on the multiband Jahn- 
Teller effects in the el-ph interaction in fuUerides com- 
pounds. We found that a band degenerate model leads to 
an almost cancellation of vertex diagrams enforcing the 
validity of ME theory, as observed by Quantum Monte 



Carlo DMFT data. On the other hand, we also show that 
realistic band structure calculations suggests that the el- 
ph interaction in A^Cqq compounds is mainly dominated 
by intraband scattering within a single band. Note that 
similar conclusions are not expected to apply for electron- 
electron (Hubbard) interaction where relevant scattering 
processes are not restricted in an energy window around 
the Fermi level. Note also that electronic and vibrational 
disorder could in principle strongly mix intra- and inter- 
band el-ph scattering. 

The effective role of the non degeneracy on the metal- 
insulator transition driven by the Hubbard repulsion and 
the inclusion of possible crystal splittings and disorder 
will be the future developments of the present work. 
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